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■ Abstract 
^ . 

_ We consider the free fall of a sphere above a wall in a viscous incompressible fluid. We 

■ investigate the influence of boundary conditions on the finite-time occurrence of contact 
between the sphere and the wall. We prove that slip boundary conditions enable to 
circumvent the "no-collision" paradox associated with no-slip boundary conditions. We 
also examine the case of mixed boundary conditions. 
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1 Introduction 



Understanding the close-range dynamics of particles inside a fluid is crucial to several fields 
of application, such as rheology, sedimentation or slurry erosion. A related problem is to 
compute the drag induced by nearby solid bodies: it is important to lubrication theory, or to 
00 ■ the control of microswimmers, among others. In the case of a viscous incompressible fluid, a 

On I standard way to model the fluid-solid interaction is through the coupling of the Navier-Stokes 

' equations for the fluid part and the conservation of momentum for the bodies. However, such 

■ model reveals unrealistic features: in particular, it can overvalue the drag force induced on 

(•~^ i bodies that are close to collision. Consider for instance the free fall of a rigid sphere over a 

I wall, and assume standard no-slip conditions at the solid surface of the sphere and the wall. 

Then, the model predicts that no collision is possible between the sphere and the wall, no 
matter their relative density and the viscosity of the fluid. This is a "no-collision paradox": it 
^ I has been recognized at first for reduced models (Stokes approximation, see [5l[6]), and more 

recently for the full constant-density incompressible Navier Stokes system \13\ [T9\ [2T] . These 
results show the need for a more accurate model of the fluid behavior in the "small-distance" 
regime. Several refinements were proposed in the past [H [HI [23, see the introduction of 
[26j for a review. In this paper, we investigate the infiuence of boundary conditions at the 
fiuid-solid interfaces. 



*Institut de Mathematiques de Jussieu, CNRS UMR 7586, Universite Paris 7, 175 rue du Chevaleret, 75013 
Paris. Mail: gerard-varet@math.jussieu.fr 

^Ceremade, CNRS UMR 7534, Univerisite Paris Daupliine, Place du Marechal de Lattre de Tassigny, 
75775 Paris Cedex 16 and Institut Camille Jordan, CNRS UMR 5208, Universite Lyon 1, 43 boulevard du 11 
novembre 1918, 69622 Villeurbane Cedex, Mail: hillairet@ceremade.dauphine.fr 

^Institut de Mathematiques de Jussieu, CNRS UMR 7586, Universite Paris 7, 175 rue du Chevaleret, 75013 
Paris. Mail: wangc@math.jussieu.fr 



1 



Boundary conditions are constitutive equations that can be derived theoretically or ex- 
perimentally. A popular general type of boundary conditions was proposed by Navier in 
the early time of analytical fluid dynamics. This Navier law states that, on the one hand, 
the normal component of the velocity is continuous at the fluid/solid interface (imperme- 
ability condition) and that, on the other hand, the amount of slip in the tangential part 
of the velocity is proportional to the tangential part of the normal stress exerted by the 
fluid on the boundary. The proportionality coefficient is called "slip length". Assuming the 
slip length vanishes {i.e., assuming no-slip boundary conditions) suffices to make theoretical 
and experimental results to coincide in numerous cases. However, developments such as the 
"no-collision paradox" lead to question this assumption [23]. We refer the reader to [2l I28j 
for a further description of experimental and theoretical validations of the Navier boundary 
conditions with slip. We propose here new insights on the contact-problem, by considering 
the effect of Navier boundary conditions. We shall study this effect in the context of the full 
Navier-Stokes system, which is, up to our knowledge, new. 

We consider a simplified case where only one homogeneous solid is moving in the fiuid. We 
denote S{t) C M.^,F{t) C the solid and the fluid domains at time t, and Q = ^S(t)UF(t) the 
total domain. We assume that the fluid is governed by the constant-density incompressible 
Navier-Stokes equations : 

pridtup + Up ■'^up) — fJ-F^up = —Vpp — ppge^, t > 0, x e F{t), (1.1) 

div up = 0, t > 0, X E F{t), (1.2) 

where up and pp stand respectively for the velocity and internal pressure, pp the density, pp 
the viscosity in the fluid. We emphasize that we consider here the influence of gravity whose 
direction is the third vector of the basis (61,62,63), i.e., the "vertical" unit vector. The solid 
motion is governed by Newton laws. Denoting xs{t) S the position of the center of mass 
of the solid at time t, Us{t) € its velocity and ujs{t) € the angular velocity, the Newton 
laws read 

d f 

ms-j-Usit) = - {2ppD{up) - ppl)nda - msge^, (1.3) 
JdS{t) 

d f 

-(Jsujsit)) = - {x-xs{t))x{2ppD{up)-ppI)nda. (1.4) 
"I JdS{t) 

Here ps and ms — P5|S'(0)| are the density and the mass of the solid respectively. Symbol n 
stands for the unit normal vector to dS{t) pointing outside the fluid domain. The 3x3 inertia 
matrix -Jsit) is defined by 

Jsit) = ps [ {\x - xsit)\'^Id -{x- xs{t)) (g){x- xs{t)))dx. 
Js{t) 

We emphasize that the solid is homogeneous, so that no angular momentum is induced by 
gravity in (|1.4p . Given {Us,u}s, xs), the solid velocity-field us at each point x E S{t) reads 

us{t, x) = Us{t) +iosx{x- xs{t)). (1.5) 

In this paper, we complete the system with Navier boundary conditions : 

{up - us) ■ n|a5(f) = 0, {up - us) X n\Qs(t) = -2/3s{D{up)n) x n, (1.6) 
UF-n\dn = 0, Up X n\dn = -2Pn{D{up)n) X n, (1.7) 
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where f3s,Pn ^ are the shp lengths at the fluid-sohd and fluid-contamer interface respec- 
tively. We remind that assuming that f5s or /3q vanishes means assuming no-slip boundary 
conditions at the corresponding interface. 

System ()l.ip - ()1.7p is completed with initial conditions : 

S{0) = So, uf{0,x) =UF,o{x), x£n\So (1.8) 
U{0) = Uo, oj{0)=coo. (1.9) 

Equation ()1.9p may also be written : 

us{0,x) = us,o{x), xeSo- (1.9) 



where usfl is computed with respect to Uo,ujq and via (jl.Sp . When Sq C 0, and uf,o 
has finite energy, namely uf,o G L'^{0,\ Sq), standard computations show that a reasonable 
solution to ()l.ip - ()1.9p should satisfy the energy estimate : 



PF\uF{t)\'^+ / Ps\us{t)\'^ 
F(t) JS(t) 



+ fiF{2l I \D{uFit))\'' + ^ I I \uFxn\'^da + — / \{uf - us) x nfda 

Ian PS Jo Jds{t) J 



(2 f [ \D{uFm' + ^ f j 

\ Jo JF{t) Pn Jo Jai 

<-/ / pFge3-UF- I I Psge3-us + l:[ Pf\uf,o\^ + 7; I Ps\us{0)\'^ ■ (1.10) 

Jo JF(t) Jo Js{t) ^ JFo ^ J So 

where Fq = \ So- When /55 = or /3f7 = 0, the corresponding boundary term vanishes. 

Roughly, with regards to the classical theory of the Navier-Stokes equation in a fixed 
domain, one can expect two kinds of solutions for (ll.ip - (ll.9p : 

• Strong solutions, locally in time: they should be unique, with possible blow-up prior to 
collision in 3D. 

• Weak solutions: they should exist at least up to collision, and be possibly non-unique. 

In the no-slip case {(3s = (3n = 0), both kinds of solutions have been built in various contexts: 
2D/3D problems, bounded/unbounded container, smooth/singular shape of the solid bodies. 
Existence of strong solutions locally in time (prior to collision) is now well-known, see [9l 
[T2t [T3t \T7\ \32\ [33] . The existence of weak solutions up to collision has been studied in 
lilElllOlllHllMlElllSOllM]. Uniqueness of weak solutions up to colhsion in 2D is estabhshed 
in the recent paper [16] , Finally , extension of weak solutions after contact is given in [29] 
for 2D problems and [11] for 3D problems. 



In the case of slip conditions, when /3s > and (3^ > 0, much less in known. It is 
partly because genuine new mathematical difficulties arise: for instance, at the level of weak 
solutions, discontinuity of the tangential velocity at the fluid/solid interface forbids global 

bound on the extended velocity-field u as defined by: 

u = Us in S{t) , u = uf in F{t). (l.H) 
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Such global bound is a key point in the treatment of the Dirichlet case. Hence, to investigate 
the contact problem for system (jl.ip - ()1.9p . we shall rely on the recent paper [15], which 
shows the existence of weak solutions up to collision. Both the definition of weak solutions 
and the existence result given in |15j will be recalled in the next section. Let us also mention 
the work [31], related to a single body without container. 

Considering a special class of weak solutions, we will show that slip conditions allow 
collision in finite time. Hence, taking into account slight slip at the solid boundaries allows 
to clear the no-collision paradox. More precisely, we shall consider configurations obeying 
the following assumptions: 

Al. The solid S{t) ~ is a ball of radius 1. 

A2. is a smooth convex domain of {xs > 0}, flat near: 

50 D {x3 = 0, l^l < 2(5} for some 6 G (0, ^). 
A3. S{t) is axisymmetric with respect to {xi = X2 = 0}. 

A4. The only possible contact between the solid and the container is at x = 0. More 
precisely, 

dist {S{t),dn \ {x3 = 0, |x| < S}) >ds>0. 

Parameters 5 and ds will be involved in the proof. Let us comment on these assumptions. 
Their only goal is to ensure that the geometry of the potential contact zone does not change 
with time. Namely, it is an aperture between a sphere and a plane, only varying with the 
distance h{t) between the south pole of S{t) and x = 0. The whole point is then to determine 
if h{t) can vanish in flnite time or not. This simplifles greatly our computations. However, 
as will be clear from our proof, our result is genuinely local in space time: it is independent 
from the global characteristics of the solid/container geometry. 

Note that A1-A2 are just assumptions on Sq and fl. If So CC J^, and if the initial data 
{up^Q, usfi) has finite energy, article |15j ensures the existence of a weak solution up to contact 
between S{t) and di} (see next section). Assumption A3 states that the axisymmetry of the 
solid is preserved with time. Actually, it will be automatically satisfied if Sq, ^ and the initial 
velocities are axisymmetric. More precisely, in this case, one can adapt the construction 
of weak solutions from |15] by including the symmetry constraint, and obtain directly an 
axisymmetric weak solution. Finally, assumption A4 guarantees that the sphere remains 
away from the top and lateral boundaries. Again, this can be ensured by a proper choice of 
the container and the initial data. As regards lateral boundaries, it is enough by assumption 
A3 that they are at distance > 1 to the vertical axis {xi = X2 = 0}. Also, there are several 
ways to avoid merging of the solid body and the top boundary. One can for instance assume 
that the solid is heavier than the fiuid {ps > Pf), and that the initial kinetic energy of the 
fluid and the solid is small enough. As the total energy of weak solutions is non-increasing, 
this means that the increase of the potential energy through time can only be small: hence, 
the solid will not rise too much, and so will remain away from the top boundary. Again, we 
insist that our method of proof could extend to more general settings. 

We can now state our main results: 
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Theorem 1.1 (Slip Case). Assume fisil^Q, > '^''T'd Ps > Pf- For any weak solution {S,u) 
(see Definition \2. 1\) satisfying AI-A4, the solid body S touches di} in finite time. 

Theorem 1.2 (Mixed Case). Let j3s or /3q vanish. Any weak solution {S,u) (see Definition 
\2.3]) satisfying AI-A4 is global. In particular, the solid body S never touches dO,. 



We recall that the full no-slip case (when (5s = Po. = 0) was already treated in |21j . 
Theorems 1 1 . 1 1 and 1 1 . 2 1 show that one possible way to obtain realistic contact, and to circumvent 
the no-collision paradox, is to introduce slip in the boundary conditions on all solid boundaries 
of the fluid domain. 

We conclude with a few hints on the proof of Theorem 11.11 (ideas for Theorem 11.21 are 
similar). The proof elaborates on the paper [T3], devoted to a simplified linear system: one 
takes there = Rij_, while the fluid outside the sphere S{t) is governed by a steady Stokes 
flow (with //i? = 1 for simplicity): 

—Aup = —VpF — ppg^s, V • Up = 0, X £ F{t). 

In this simplified setting, the fluid and solid domains F(t) = F^-^ and S{t) = Sh(t) are 
characterized by the distance h{t) beween the sphere and the plane wall. Then, it is easily 
seen that upit) = h'{t)wh(^t) where vuh is the solution of a normalized Stokes equation, set in 
the domain "frozen" at distance h: 

-Awh + yph = 0, V • = 0, X £ Fh, [wh - 63) • n|a5^ =0, Wh- n\dQ. = 0, 
{wh - 63) X n\ds^ = -2(3s{D{wh)n) x n, 
Wh X n\QQ, = -2/3n{D{wh)n) x n. 

(1.12) 

We refer to [13] for all details. Eventually, the dynamics reduces to an ODE of the type 

h"{t) = -h'{t)V{h{t)) + PS~P^ g 

PS 

where the drag term T){h) is the energy of w^: 

V{h) = 8h{wh), £h{u) := [ |Vn|2+(-l + l) / \(u - es) x + ^ [ \u x n\\ 

JFu \Ps J JdSh Jan 

For this reduced model, the contact problem resumes to the computation of the drag T>h in 
the limit /i — )• 0. The difficulty is that there is no simple formula for the solution Wh of the 
Stokes system in Fh. 



A general approach to this problem was recently proposed in |14j . It works for no-slip 
conditions, slip ones, and for more general geometries. The starting idea is to use a variational 
characterization of the drag: one can identify I'(/i) as the minimizer of an energy £h, over an 
admissible set of fields Ah- 

V{h) = inf £h{u) = £h{wh). 
u&Ah 

Then, instead of computing the true minimizer, the point is to find an appropriate relaxed 
minimization problem, for which the minimizer can be easily computed, and such that the 
corresponding minimum is close to the exact one asymptotically in h. The choice of such 
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relaxed problem is developed in various situations in [T3] . Roughly, in the case of slip condi- 
tions, one obtains that c| ln(/i)| < D{h) < C\ ln{h)\, in agreement with older computations of 
Hocking |23j . Hence, the drag is weak enough to allow for collision, within this linear Stokes 
approximation. See also [20] for numerical insight. 



In the present paper, we manage to use ideas of the linear study, in the context of the 
nonlinear Navier-Stokes flow. At a formal level, the first idea is to multiply the momentum 
equation by w{t,x) = Wi^(^i^{x), Wh solution of (|1.12p (where this time is the container). 
After formal manipulations, one obtains 

h"{t) = -h'V{h) + g + remainder (1.13) 

PS 

where the remainder term involves the test function w and the convective derivative of u. 
However, the Stokes solution is not known explicitly. Hence, the main point in our proof 
is to use another test function, based on the relaxed minimizer introduced in \14^ . Then, one 
must show that for such test function, the remainder term is indeed small, that is controlled 
by the gravity when the solid is heavier than the fluid. Let us stress that several difficulties 
arise in this process. For instance, the relaxed minimizer introduced in the linear study is 
not accurate enough, so that we must improve it. Moreover, due to the slip conditions, the 
natural quantity is the symmetric gradient of the weak solution, not the full gradient. This 
makes the derivation of the various bounds more difficult than in the no-slip case. 

The outline of the paper is as follows. In the next section, we introduce a suitable notion 
of weak solutions for system (jl.ip - (jl.9p . in the slip case and in the mixed case. We recall the 
existence theorems up to collision available for such weak solutions, in the slip case and in 
the mixed case. The two last sections are devoted to the proofs of the two theorems. 

2 On the definition of weak solutions 

We introduce here a definition of weak solutions for (jl.ip - ()1.8p in two cases: slip boundary 
conditions at all boundaries, mixed slip/no-slip boundary conditions. We mostly follow |15] . 
Given Q a Lipschitz domain, we set: 

v^{n) ^ {^ec^{n), divv? = o}, v^m ^ win, ^ ev^iw")}, 

Ll{n) = the closure of V^{n) in L'^{n), Hl{n) = H^{n) n Ll{n), 

= the closure of V^{Ti) in H'^{n), 

n = {^ps, ips = V + (jj X X, for some yGM^,wGM^}, 



2.1 The case of slip boundary conditions: > and > 0. 

We consider here that the slip lengths f3s and (Sq do not vanish i.e., that the fluid slips on 
both boundaries. The corresponding Cauchy problem is studied in [15], where the following 
definition of weak solution is given: 

Definition 2.1. Let 17 and So Q, be two Lipschitz bounded domains of M^. Let uf,o G 
L'^{i^),usfl € TZ such that upfi • n = usfl • n on OSq. Let T > 0, or T = oo. 
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A weak solution of \1. 1]) - I[T79\) on [0, T) is a pair {S, u) satisfying 

1) S{t) is a bounded domain ofR^ for all t £ [0,T) such that 

Xs{t,x)^lsit){x)eL^{{0,T)xn). 

2) u belongs to the space 

ST = {ue L°°(0,r;L2(fi)), there exists up S Lf^^{[{),T)- Hl{n)), us G L°°(0,T;7^) 

such that u(t, •) = upit, •) on F(t), u{t, •) = us{t, •) on S{t), for all t £ [0, T)} 



where F{t) = n\S{t). 

3) For all ip belonging to the space 

TT = {ip£ C{[0,T];Ll{fl)), there exists ifp S V{[0,T);V^{Ti)),ips G P([0,^);7^) 

such that ip{t, •) = (prit, •) on F{t), ip{t, ■) = ^ps{t-, •) on S{t), for all t G [0, T]} 

there holds 

/ PFUF ■ dtifF - / PSUS ■ dt(ps - / PFUF ® UF ■ VifF 

JF{t) Jo Js{t) Jo J F{t) 

+^ / {uF xn) ■ {ifF xn) + ^ / {{uF - us) x n) • {{(pF - (ps) x n) 



Pn Jo Jan Jo Jds(t) 

(■T r f-T r rT 



+ / 2iifD{uf) ■■ D{(Pf) = - / pFge3(pF - / Psges^fs 

Jo JF{t) Jo JF(t) Jo JS{t) 

+ / Pfuf,oV>f{0,-) + Psus,o^s{0,-), (2.14) 

J Fa JSu 



'Fo J So 

4) For allije V{[0,T);V{n)), we have 

Jo Js(t) Jo Js{t) JSo 

5) We have the following energy inequality for a. a. t £ (0, T) ; 



1 

2 

h PF 



/ pF\uF{t)\'^+ / ps\us{t)\ 
JF(t) Js(t) 

{2 [ f \D{uF{t))\^ + ^ f f \uFxn\^da + ^ [ [ \{uf - us) x n\^da] 
\ Jo JF(t) Pn Jo Jdn Ps Jo Jds{t) J 

- pFge3-UF- / Psges-us + l; PF\uF,of + 7; Ps\us{0)\'^ ■ 

Jo JF{t) Jo Js(t) ^ Jf{o) ^ Js{o) 

(2.15) 



We refer to [15j for a detailed discussion of this definition. Note that it is restricted to 
solutions up to collision, through the condition S{t) d fi. We recall that item 3) is the identity 
one gets by multiplying (jl.ip by a divergence- free test-function ipF and then combining with 
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(jl.3|) - (|1.4|) . Item is the weak form of the transport equation satisfied by xs^ the indicator 
function of S{t) as defined in 1). Finally, 2) summarizes the regularity that is implied by the 
energy estimate ()2.15p . Note that up and us share the same normal component, so that the 
]^ooj2 ]-,Q^]2d on Up and us yields a global L°°L'^{yt) on the extended velocity field u (see 
(jl.lip for the definition of u). Moreover, the L°°L?' bound on up and L^L^ bound on D{up) 
imply an Lf^^H^ bound, through Korn inequality (that is valid as long as the solid is away 
from the boundary of the container). 

Prom [15], there is a weak solution for the system up to a contact between S{t) and dO,: 

Theorem 2.2. Let Q and 5o (s ^2 two C^'^ bounded domain ofW^. Given (/35,/3q) G (0, oo)^, 
and initial data up^ G L'^{^l),us,o G TZ satisfying: 

upfl ■ n = usfi • n, on OSq, 

there exists T G R!j_U{oo} and a weak solution of U.1\ )- [T7^) on [0, T). Moreover, such weak 
solution exists up to collision, that is either we can take T = oo, or we can take T G R!j_ in 
such a way that 

S(t)'£n, for all iG[0,T), and Urn dist(dS(t), dQ) = 0. 

t-fT- 

The aim of Theorem 11.11 is to prove that T < oo under assumption A1-A4. 
2.2 The mixed case: (3s = and (3^ > 0. 

We proceed with the case where slip is imposed on one of the boundaries only. For simplicity, 
we consider that no-slip boundary conditions are imposed on dS{t) i.e. I3s = and I3q > 0. 
To our knowledge, this particular system has not been treated explicitly in previous studies. 
Nevertheless, as slip is allowed on the exterior boundary of F{t) only, it is straightforward to 
adapt the method of [10] or |25j . for instance, to tackle a Cauchy theory for weak solutions 
before contact. In particular, the extension trick (jl.lip yields a sufficiently smooth velocity- 
field in this case so that we can reduce the full system to a global weak formulation on this 
extended velocity field. This reads as follows: 

Definition 2.3. Let and Sq C 0, be two Lipschitz bounded domains o/M^. Let up^ G 
L'^{Q,),usfi G TZ such that up^ ■ n = usfi ■ n on OSq. A weak solution of iLl\) - [T79\) is a pair 
(5, u) satisfying 

1) S{t) 1^ Q. is a bounded domain o/M'^ for all t G [0,T) such that 

Xs(t,x)^l5(i)(x)GL-((0,r)xf]). 

2) u£ L°°(0, T; (fi)) n L'^{0, T; H^in)) satisfies 

w(t,-)ls(,3 for a.a. t e {0,T) . 

3) For all belonging to 

TP = We V{% T) X 0), (^(t, Oi.^,) G n for all t G [0, T)} 
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there holds 



p{udtip + u®u: D{ip)) + I I 2fiFD{u) : D{ip) 

Jo Jn 

pgesip- / pQUQipiO,-), (2.16) 
Ju 



/3n 



+ I I {u X n) ■ {if X n) 

/o Jdn Jo Jn 



where 



p{t, x) = PS ls{t) {x) + pplpit) (x), for (t, x) G (0, T) x n, 
Po{x) = pslsoix) + PF^Foix), forxeQ, 
Uo{x) = us,o '^Soix) + upfi'i^Foix), for X £n, 

4) For alli^e V{[0,T);V{Q)), we have 

- r [ dt^- r [ u-vij= f v(o,-), 

JO JS(t) Jo JS(t) JSo 

5) We have the following energy inequality for a. a. t G (0,T) ; 

"t f ^ ft 



\ [ p\n{t)\' + 2pF [ [ \D{u)\' + l- 
^ Jn Jo Jn 



\u X nPdcr 



Jdn 
< - 



f I pge^u + l I /9(0,-)|^o|'. (2.17) 

Jo JF{t) ^ Jn 



Following the method of proof of |T0], one obtains: 

Theorem 2.4. Let and S'q (s two C"^'^ hounded domain o/M'^. Given (3q £ (0,oo) and 
initial data uf,o G L'^{^),usfl G TZ satisfying 

uf,o ■ ^ = ""5,0 ■ on dSo, 
the same conclusion as in Theorem fHH holds. 



3 Proof of collision with slip conditions 

This section is devoted to the proof of Theorem ll.li We assume fis and /3q are fixed non- 
negative slip lengths, and ps > PF- Let (S, n) be a weak solution over (0,T), as given by 
Theorem 12.21 We consider a solution up to its " maximal time of existence" : that is T = oo 
if there is no collision, or T < co the time of collision. We assume that this weak solution 
satisfies A1-A4 for t G [0,T). We denote hy h = h{t) the distance between the south pole of 
S{t) and X = 0. By assumption A4, when h{t) is small enough, it is equal to the distance 
between S{t) and the boundary d^l. Our goal is to prove that h{t) goes to zero in finite time. 
In other words, T < oo, that is collision occurs in finite time. 

To build up our proof, we first need some general considerations on the connection between 
the weak formulation (j2.14p and the distance h, for an arbitrary weak solution satisfying A4. 
We then design a special test-function in order to prove contact in finite time. 
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3.1 Some features of the weak formulation fl2.14p 

Let (5, u) a weak solution defined on (0, T) and 99 G 7^ of the form 

f = C{t)^hit) with ips = C{t)e3 . (3.18) 

for some functions ( G P([0, T)) and h 1— )• ifh to be made precise. Integrating by parts, we 
get that for any q G V{[0,T) x O), for a.a. t G (0,r) 



2 / D{^f) ■■ D{uf) 

JF{t) 

= [ {2D{^f) - ql) ■■ 
JF{t) 

= — {AifF — Vq) ■ uf + / {2D{ipf) — ql)n ■ UFda 

JF(t) JdF(t) 

= — I {/\(pF — '^q) ■ uf + / {2D{ipf) — ql)n ■ UFda + / {2D{ipf) — ql)n ■ UFda 

JF{t) Jdn JdS{t) 

= — [AifF — ^q) ■ uf + / {2D{ipf) — ql)n ■ UFda 

+ / {2D{ipF) - ql)n ■ {uf - us)da + / {2D{lpf) - ql)n ■ usda, 
JdS{t) JdS{t) 

Chosing also q of the form q = C(*)9h(t); the last term can be calculated by 

/ {2D{ipF) - ql)n ■ ugda = /i'(i)C(i) / (21?(<^h{t)) " qh{t)I)n ■ esda ^ Q{t)h' {t)n{h{t)) . 
Jas{t) Jds{t) 

We emphasize that the domain S{t) is completely fixed by the value of h{t) (we shall denote 
Sh{t) in what follows) so that the last integral does depend on h{t) only. Combining all the 
equations with (j2.14p . we find 

^lF I C{t)h'{t)n{h{t))dt 



Jo 

/ PFUF ■ dtifF + / / PSUS ■ dtfS + / PFUF UF : VipF 

JF{t) Jo Js{t) Jo JF{t) 

- PFgesfF- / / Psge3^s+ / PFUF,oVFiOr) + / Psus,o^siO,-) 

Jo Jfu) Jo Js(t) Jfo Jso 



^t) Jo JS{t) JFo JSo 

^ {uF >^ n) ■ {(fF X n) - ^ / {{uF - us) x n) • {{ipF - fs) x n) 

Pn Jo Jp PS Jo Jasit) 



+ PF 



( [ [ {ALpF-Vq)-UF- [ I {2D{(Pf) - ql)n ■ UFda 

\ Jo J F(t) Jo J do. 

[ [ {2D{^f) - ql)n ■ {uF - us)da] (3.19) 
Jo Jas{t) J 
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This identity enables to describe the time evolution of h through the function n that we have 
to choose now. To compute n{h), we shall take advantage of the following identity: 

n{h) = / {2D{iph) - qhl)n ■ (fhda + / {D{(ph) - qhl)n ■ (es - <fh)da 
JdSh JdSh 

= / {^Vh - ^Qh) ■ Vhdx + 2 D{iph) : D{iph)dx 

{2D{(fh) - qhl)n ■ (fihda + / {D{(ph) - qhl)n ■ (es - '^h)da. (3.20) 
an JaSh 

where = Q.\Sh- This identity holds true for all h such that Sh CC Q.. 
3.2 Construction of test-functions 

We construct now the test function and the pressure associated with the solid particle 
Sh frozen at distance h. It will be defined for h G (0, with Hm — ^'^'9t&[Q,T) h{t). When 
/i — 7- 0, a cusp arises m. = Vt \ Sh. This cusp is contained in the domain : 



^h,r ■= {{xi,X2,X3) e Fh, X3 < - + h and |(xi,a;2)| < r} , 

for arbitrary r € (0, ^). In particular, in all forthcoming estimates, we shall pay special 
attention to the region 0,h^2S, where 6 was introduced in assumption A2. 

To compute test-functions and associated pressures, we introduce cylindrical coordinates 
{r,9,z) associated with cartesian coordinates (xi,X2,X3) in such a way that X3 = z. The 
local basis reads {er,ee,ez)- First, we introduce an approximate solution of the steady Stokes 
equation in the domain il,h,25- It will satisfy locally the slip boundary conditions (jl.6p - (|1.7p . 
with Us replaced by e^: 



1 

where 
Here 



-dz'^ r, 



h + 7s (r) 



xi, -dz^ [ r, . ) X2, -dr (r, 



sir) = 1 - VT 

Pi{r) 
P2{r) 



^{r,t)=Pi{r)t + P2{r)t^ + Pzt^ 



6(2 + 05) 



12 + A{as + ap) + agap 

3(2 + as)ap 
12 + 4(05- + ap) + asap ' 

2{as + asap + ap) 
12 + 4(05 + ap) + asap ' 



where 



ap 



h + jsjr) 



as 



^ + 2)(/i + 7,(r)). 



(3.21) 

(3.22) 
(3.23) 
(3.24) 

(3.25) 
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This approximation is inspired by the computations in [T3]. We emphasize that the value 
oi as does not match exactly the one in [T3]. This modification is required so that further 
boundary estimates hold true (see p.39p ). We stress that notation d^^{r, ), resp. 

). With this 



*fr, z) = <&fr, 



dr'^ir, h^^), stands for 5^'*(r, z), resp. 5;'^'(r, z), 

structure, iph is divergence-free. Moreover, the function is constructed so that it satisfies 
the following boundary conditions: 

• for z = 



^>(r,0) = 0, 



h + 7s (r) 



h + 7s (r) 



for z = /i + 75(r) 

$(r,l) = 1 



h + 7s (r) 



+ ( + 2 ) 5,<5 ( r, 



/35 



h + 7s (r) 



In the local cylindrical basis, the approximate solution (ph reads: 



-dz^ r. 



rcr + -dr { ^ (r, 



h + 7s (r) 



(3.26) 
(3.27) 

(3.28) 
(3.29) 

(3.30) 



It 



/i + 7s(r)y r 

In particular, if we remark that n satisfies n = —rcr + {1 + h — z)ez on dSh H ^h,i/2 
n = 63 on di} n f^h,i/2) obtain (see the introduction of Appendix El for further details) : 

i^h - 63) ■n\ds^ndnh,2s = 0' 

• ^^|onne^^,2^ = 0' (3-31) 
^fe X n|aQne!^h,2i = -2/3n(-D(^h)n) X n|af7n9n,,,2i- 
Thus, up to its extension to the whole of 0, the function ipf^ is a good candidate for applying 
the computations of the previous section. Such extension iph should satisfy (ph\sh = ^- 
remains to define it into the whole of F^. 

To this end, we introduce a smooth function <j) = <j){x) with compact support satisfying 
0=1 on a ^-neighborhood of the unit ball -6(0, 1) , 
(j) = outside a d^-neighborhood of B{0, 1). 
We recall that ds was introduced in A4. We also define a cut-off function x 

X = 1 in (-6, 6f, X = outside {-26, 25 f . 
With these cut-off functions, we set (in cartesian coordinates): 



[0, 1] s.t. 
(3.32) 



-X2dz 



(1 - x{xm - (1 + h)e3) + x{x)Hr, j^^) 
(1 - x{x))^ (x - (1 + h)e3) + x{x)Hr, j^^) 



, -rdr 

63 in Sh- 



r2((l - xix))^ ix-il + h)e3) + xix)'^ir, 



)) 



in Fh, 



(3.33) 



By the definition of (ph, we get 
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Lemma 3.1. The map {h,x) i— )■ ^h{^) is smooth on all domains: 

{{h,x),h£ (0,/iA/),x G Fh}, {{h,x) G {0,hM),x G Sh} , {{h,x),h€ %hM),x£ F^X^h^s}- 
Furthermore, there holds: 



and 



iph ■ n = on 



fPhlFh ■ = €3 ■ n ondSh 



(3.34) 



(3.35) 



This comes from direct calculations, that we leave to the reader for brevity. We now state 
some refined estimates, to be used in the proof of Theorem II. 1[ The whole point is to have 
an accurate bound on the blow-up of the test function as h goes to zero. 



Lemma 3.2. There exist constants < c < C so that, when h < hM ■' 

\\'^h\\L2{Fh) ^ 
W^^hWh^F,) < C\lnh\, 
\\D{iph)\\l2^p^^ > c|ln(/i)|, 

and 

\\2l3s{D{iph)n) xn+{iph- eg) x n\\L2(^Qs^) < C, 
Furthermore, fori = 1,2,3, denoting by exponent i cartesian coordinates, there holds: 

dhVhdA\l^9Qndn^_,) ^ C\lnh\, 



\\{h + ^s{rm 



dh^\{xi,X2,s)ds\\L2^^^^^) < C, 



(3.36) 
(3.37) 
(3.38) 

(3.39) 

(3.40) 
(3.41) 
(3.42) 



and for {i,j) G {1, 2, 3}^ and x G 



/ {di(p{{xi,X2,s) + dj(pl{xi,X2,s))dt 

J z 



h+7s(r) 



{diipl{xi,X2,s) + djipl{xi,X2,s))dt 



< C, 



< C. 



(3.43) 
(3.44) 



The proof of this lemma is postponed to Appendix Rl 

When computing (j3.19p . we introduced a pressure q. The aim was to smoothen singular- 
ities in Aifh which must arise when /i — >• 0. The construction of this pressure is the content 
of the next lemma: 

Lemma 3.3. There exists a map (h, x) i— t- qh which is smooth on the domains 

{{h,x),h G (0,/iAf),x G Fh} and {{h,x),h G [0,/iA.f),x £ Fh\ ^h,5} 
and which satisfies : 
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there exist constants < c < C < oo such that : 



c\lnh\< (2D{iph) - qhl)n ■ {es - iph)da - {2D{iph) - qhl)n ■ (phda < C\\nh\; 
JaSh J on 

• there exists a positive constant C such that, for any h G (0, /ijv/) and v G H^{Ffi) 
satisfying v ■ n = on d^l, we have: 

[ {^^h-^qh)v <C{\\D{v)\\L2^p^) + \\v\\L2i^Qn))- (3.45) 

The proof of this lemma is also postponed to Appendix 1X1 
3.3 Conclusion of the proof of Theorem 11.11 

We are now in a position to prove Theorem II. 1[ Combining the energy inequality (j2.15p with 
the classical identity 

- /* / ppge. -up- f f psge. ■ us = ^iP^—P^{h{t) - h{0)) 

Jo JF(t) Jo Jsit) -3 

we obtain easily (remember that ps > Pf)- 

sup / \uit)\'+ r [ iDiupm^ 
{o,T)Jn Jo J F(t) 

+ I I |nFXnpdcr+ j [ {{up - us) ><■ n\'^da < Cq. (3.46) 
Jo Jan Jo Jds{t) 

This energy bound will be of constant use. 

As mentioned before, we take a special test function in the variational formulation (j2.14p . 
Namely, we take (p{t,x) = C(t) (^/,(t)(3;), with the field iph built in paragraph l3.2[ From Lemma 
13.11 we have that G 7t so that identity ()3.19p holds true. The whole point is then to bound 
properly the r.h.s. of (j3.19p : 



RHS =11 ppUF ■ dtipF + / / Psus ■ dtfs + / Pfuf ® up : V(pp 

Jo JF(t) Jo J S(t) Jo JF{t) 

- I \ I PFgez^F+ / Ps9ezVs\ + \ / PFUF,o^F{^r) + / Psus,o^s{^r)] 

Jo \JF{t) Js{t) J \JFo JSo / 

rT 



pp [ [ {up X n) • ( -^{<fF X n) + {2D{lpp) - ql)n] da 

Jo Jan \m J 

pF / {{uF - us) X n) ■ (-^{<fF - fs) xn) + {2D{(fF) - ql)n] da 

Jo Jasit) \PS J 



+ PF / {A(fF - Vq) ■ up 
Jo JfU) 



'Fit) 

h+h + h + h + h- fiFih + I7 + h), 



14 



where we have apphed the continuity of normal traces on dS{t) and in the boundary 
terms. 

• We first deal with 

rT r r-T 



3 



Jo Jfu) Jo 



We introduce the functions x i— )• Pz^p^x) for (i, j) S {1,2,3}^ defined by: 

{dnpl + dj(pl){xi,X2,s)ds , Vr<l, z < h + 'j^ir), 

Pz^h^ix) = / {diipl^ + djipl^){xi,X2,s)ds otherwise. 

Let us stress that PzV^^if is smooth in Ff^, notably across r = 1: it follows from the cancellation 
of the symmetric gradient dup'j^ + dj(pl^ near (r = l,z = 1 + h) that both expressions coincide 
there. 

Moreover, by (j3.43p . there exists a constant C independent of t such that, for all G 
{1,2,3}2, there holds: 



\PzVh\xi,X2,X3)\ + \{h + 'ys{r))diPz(pl\xi,X2,X3)\ < C, yx£nh{t),s, (3.47) 



\Pz^]f{xi,X2,X3)\ + \d,Pzy^]f{xi,X2,X3)\ < C, Vx G F{t)\nhit^^s- (3.48) 



Then, by performing integration by parts, we obtain (summation on indices i and j is im- 
plicit): 



hit) = - [ u^FdziPz^l'] = f dsiiJ^FU^Pz^if - [ u^FU^Pzipifn^da 
JF{t) J F(t) Jan 

= [ {{dsv'F + diU^^F + id3U^ + d,u''F)u'F)Pz^h' 
JF[t) 

- / {diU%u^F + djU%u'F)Pz'pl^ - / (u^n^^)P^(^*^^n3dc7 
JF(t) Jan 

{{dsu^ + diul)u^ + (93< + djul)v^F)Pzvif 

) 

+ / ({diui, + djU^F)uFPzP>if + u%u^FdiPz<pl^ + ulu'pdjPz^'f 
Jfu) ^ ' 



F{t) 

{upu^n^ + u^Ufh' + UFUFn^)Pz(p]fda. 



an 



In this last identity, we apply (|3.47p - (j3.48p to obtain : 



{upU^F'''^^ + Up'^F'^'' + UFUFn^)Pz(p'ifda 
an 



<c\\uFmh 



(an) 
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Concerning the other terms, we have for instance : 



/ 

JFit) 



< 



F{t)nf7,(,)., h + 7s(r) 
+ 



< 



h + 'js{r) \\^F\\LHF{t)) \\{h + -fs{r))diP^ipl^ llL°°(F(f)nn„(t),,)) 

HuU\L^iF{t)n{\{xi,X2)\<{l+ds)})hF\\L^(F{t))\\diPz'fl^^ 

As Up{x) = for X G 2:2)1 < 1 + ^5,^3 = 0}, a straightforward Poincare inequahty 

yields that there exists a constant C > independent of time so that: 

\WF\\L'2{Fit)n{\(xi,X2)\<{l+ds)}) < (||93^^FllL2(F(t)) + \Wf X ""IlLacan)) • 
Introducing ([3^ - (j3^ then yields : 

1/3(01 <C ||-D(uF)||^2(p.(t)) + P(MF(t))||L2(F(t))II^^F(i)||L2(F(t)) + hF(i)|li2(gt^) , 

SO that (|3.46p implies there exists a constant C depending only on initial data for which 

I/3I <C(1 + Vt)||C|Uoo(o,t). 

• To compute Ii, we split: 

h= / PFUF-dtfF = / / pFUFC'{t)fh{t)+ / PFUFC{t)h'{t)dh(pF = Il+Il- 

Jo JF{t) Jo JF{t) Jo J F{t) 

Applying (|3.36|) yields : 

l^ll < C'||C'llLi(0,T)ll^llL°°(0,T;L2(n)) < Co || C' II Li (0,T) • 

Then, we have = pp Jq C{t)h'{t)ll{t) where, for aU t G (0,r): 



7^1 Jm 



To bound this term, we define as previously Pz'^\ which satisfies, for all i = 1, 2, 3 : 

r-/i+7s{r) 



Pz'fhi^) = J dhfl{xi,X2,s)ds , Vr<l, z<h + js{r), 

Pz'^\{x) = I dh^\{xi,X2, s)ds otherwise. 
Again, Pz'-p\ continuous in F^, notably near r = 1, as dh'-Ph vanishes near (r = 1, z = l + 
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With (|3.4ip together with (|3.42p . there is a constant C such that: 

\\Pz<fh\\L^{F{t)) + + ls{r))diPz<f\\\L2^F{t)nQn(t),2s) + \\diPzVh\\L\F{t)\n^^,)^^s) ^ (3-49) 

\\PML^ian)<C\Hh)\'/'- (3.50) 

From now on, we restrict impUcitly the integral /^(t) to |(xi,X2)| < (l + ds) as dh^fh vanishes 
elsewhere. Performing integration by parts (summation over index i is implicit), we obtain : 



F{t) 



F{t) JdQ 

< / id3U^ + diul)P,ipi+ / uldiP.ipl- / (u^pn' + ii^Fn^)P,ipi 

JF{t) JF(t) JdD, 



Relying on (|3.49p - (|3.50p and (|3.46p . we might then reproduce the arguments we used to 
compute I3, resulting in: 



\4\ < c\\a 



\h' {t)\\\iih\^\\up X n\\L2^Q^) + VT^ . 



At this point, we remark that, given the computations of [141 Section 3.2], | \n{h{t))\^/'^ is the 
minimum dissipation energy of functions v G H^{F{t)) satisfying v ■ n = e-^ ■ n on dS{t) and 
u • n = on do,. Consequently, there exists a constant C independent of time such that: 



ln(/i)|2<C / \D{up{t))\'^ + I \upxn\'^da+ \{up - us) x n\'^da 

\JF{t) Jan Jds{t) , 



and (|3.46p implies: 



\h' {t)\\ In h\ 2 \\uF X n||i2(0Q) < C. 



Next, we recall that Jp/^) pFgG3(Ph{, 



Anppg 



which implies: 



PFges^F - I I Psges^s = ^""^^^ ^'^^ C(t)dt. (3.51) 

JF{t) Jo Js{t) "J Jo 

Then, we apply ()3.3ip together with (j3.46p . It leads to 



C(t) / {up X n) ■ ( ^{iph X n) + {2D{iph) - ql)n ] da 
Jan \Pn 

C{t) I (upxn)- (^{iphXn) + 2D{tfh)n]da 
10 Jan \Pn 

< c[ m\hFxnh2^an)<C'\\aL^^o,T)VT, 
Jo 
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• By (pIHU]) together with (ICTjl : 

I/7I = [ [ [(uF-us)xn)- (^{Lph-e3)^n) + 2D{iph)n]da 
Jo JdS{t) \PS J 

< C [ |C(t)|||(nF-'U5) xn||i2(g5(j))||(¥?,,-e3) xn) + 2/3s-D(95h)?^||L2(e5(i))dt 

Jo 

< C'\\C\\l^^o,t)Vt. 

• As for Is, applying (j3.45p . we have that 

Is= [ [ {AipF-Vq)uF < c[ \Cit)\{\\D{uF)\\L2^F{t)) + \WF\\LHdn))dt 
Jo JF{t) JO 

< C'||C||loo(o,t)^^. 

• EventuaUy, we get easily the following controls: 

\h\ < IIC'IIli(o,t)) h<C. 
Combining all these estimates to bound RHS, we end up with: 

^iF [ h'{t)C{t)n{h{t))dt 







iiripF - ps)g 



r C(t)dt + 0(C7o(||CIU-(o.T) + IIC'IIli{o,t))(1 + Vt)) (3.52) 

JO 



3 

We conclude by introducing a sequence Cn in ^([0,?")) such that 

Cn l[o,T) pointwise, with ||Cn||L°° <C, ||CIIli{o,t) < C. 

A possible choice is C,n{i) = l[o,T)(i) — x{f^{T — t)), for some smooth x compactly supported 
with X = 1 near 0. We take C = Cn in p.52p . and let n go to infinity, to get 



PF 



n{s)ds + ^AEI^IEIIt < Co(l + Vt). (3.53) 
Jh{o) 3 



Introducing the bounds of Lemma 13.21 and Lemma 13.31 into p.20p . we show there exists a 
constant C such that: 

n{h)<C\\nh\ yh<hM, 

which implies that there exists a constant C for which: 

h 

n{s)ds > —C , V /t, < Hm ■ 

HO) 

Hence, we obtain a contradiction in ()3.53p if T is large enough compared to C and Co which 
are fixed by /ia/ and the initial data only. 
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4 Proof of Theorem 11.21 



The proof is very close to the one performed above for shp boundary conditions. Therefore, 
we only sketch the proof, insisting on its few specific features. 

Computation of the drag. Let (5, u) be a weak solution given by Theorem 12.41 on [0, T). 

As before, we consider a solution up to its "maximal time of existence": that is T = oo if 
there is no collision, or T < oo the time of collision. We assume that {S, u) satisfies A1-A4 
over [0, T). This time, we want to show that there is no collision, which means that h does 
not go to zero in finite time (and thus T = oo). We shall consider a test function ip G of 
the form 

f = C{t)fh{t) with = C{t)e, . (4.54) 

associated with a pressure q € T>{[0,T) x Q), of the form 

Q = C{t)qh{t)- 

The values of (ph and qt are given below. As in the previous section, we shall use functions 
C = Cn approximating 1(o,t)- Here and in what follows, we keep the convention that = 
^plF(t) aiid ifs = f^s{t) ■ With computations similar to those of the previous section, we 
obtain: 

HF [ C{t)h' {t)n{h{t))dt 



Jo 

/ pfuf • dtfF + / / pfuf UF : V(pF 

JF{t) Jo JF{t) 

l-T f f-T 



/ / PFgez^F - / Psgdz^S + / PFUFfi^F^O, •) + / PSUSfi'^siO, 

Jo J F{t) Jo Js{t) JFo JSo 

T 

(uf X n) ■ {ifF X n) 

p 



PF 

/5n Jo 



( [ [ {/^LpF-Vq)-UF- f [ {2D{ipF) - ql)n ■ UFda] (4.55) 
\Jo Jpit) Jo Jan J 



+ PF 

Again, we compute n{h) applying the identity: 
n(/i) 



/ - Vqh) ■ iPhdx + 2 D{iph) : D{^ph)dx 

Jfh Jfu 

- [ i2D{iph)-qhI)n-iphda. (4.56) 
Jdn 

Construction of the test-function. For the function (x, h) i— t- iph{x), we keep the struc- 
ture of (|3.33p with a test-function ipfi still of the form (|3.30p . The only point that differs is 
the choice of the coefficients in the polynomial We set now: 

^i(^) = IT ^2(r) = — P-sir) = -2— , 

4 + ap 4 + Qp 4-|-ap 
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with ap given by (j3.25p . We remark that this is equivalent to taking the hmit 05 — )• oo in 
our choice of test-function in the previous section. With this choice, the function <I> satisfies 
boundary conditions: 



for z = 



$(r,0) = 0, (4.57) 
77-^1 - ^^-^ f^' hJ(^ ) = (4.58) 



for 2; = /i + 75(r) 



$(r, 1) = 1 (4.59) 

d,^{r,l)=0. (4.60) 

Similarly to the previous case, we provide also the following technical lemmas : 
Lemma 4.1. The map {h,x) 1— t- iph{x) is smooth on all domains : 

{{h,x),h G {0,hM),x £ Fh} , {(/i,x) G {0,hp_[),x G Sh} , {{h,x),h G [0,hM),x G Fh\^lh,s}- 
Furthermore, there holds: 

iph-n = on on (4.61) 

and 

^h\Fh = es on dSh- (4.62) 
Lemma 4.2. There exist constants < c < C so that, when h < Hm ■' 

MmF,) < C, (4.63) 

c 



\\yvh\\i2^F,) < J, (4-64) 

\\D{iph)\\l2^F^) > (4.65) 
Furthermore, fori = 1,2,3, denoting by exponent i cartesian coordinates, there holds: 

l'h+-fs{r) 

II / dh(pldz\\L2(^onndni,s) < C\lnh\, (4.66) 
Jo 

dh^ldsh2^n,,s) < C, (4.67) 

J z 

\\{h + -fs{r))di J dhVp{xi,X2,s)ds\\L2(^Q^ ..) < C, (4.68) 

and 

• for i G {1, 2, 3}^, and x G ^h,5- 

Cr 

di^l + dsifl < (4.69) 
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for £ {1,2}2, and x £ n^^g: 



ft+7s(r) 



{dnpl{xi,X2,s) + djipl{xi,X2,s))di 



{h + js{r))di I {di(pl{xi,X2,s) + djipl{xi,X2,s))ds 

2 



< C, (4.70) 

< C. (4.71) 



Lemma 4.3. There exists a map {h,x) i— )■ qh which is smooth on the domain {{h,x),h € 
{0,hM),x G Fh} and which satisfies : 

• there exists a constant positive constant C such that : 



an 



{2D{iph) - qhl)n ■ (fhda 



< C\lnh\; 



there exists a positive constant C such that, for any h G (0, /im) and v G H^{Q) 
satisfying u G TZ and v ■ n = on dil., we have: 



[ iAiph-Vqh)v <C{\\Div)\\L2^F^) + \\v\\L2^n))- 
More details concerning these lemmas are provided in Appendix iBl 



(4.72) 



Conclusion. Introducing the results of Lemma 14.21 and Lemma l4.3t we compute again the 
right-hand side of (|4.55p with the help of energy estimates satisfied by u. We only point out 
two differences. First, there exists a universal constant C so that 

\\'^F\\L2{dS{t)) = \\us\\L2{dS{t)) < C\\us\\L2(S{t)) < C\\u\\l2(^^), 

by taking into account that us is a rigid vector field. Second, one term requires a new 
treatment (we drop summation over index i G {1,2,3} for simplicity): 



/ 

JF{t) 



U'puy{diip% + d2,ipF) 



u'i 



X- 



-u'p{ls{r) + h){d,^l + d^^'p) + 0{\\uF\\l2) 



'Fm^H(t),2, ls{r) + h 

< \\u'p{-fs{r) + h){di^% + 93V^F)llL^(Fft)nn,(,),2,)ll ^^^^^^^^ llL2(Fft)nn„(t).2.) + OiWupWl^) 

< C\\D{up)\\p2(^p^^ty^\\xu'p{-is{r) + h){diipl + 93VF)llL2{F(t)nf7h(t),2i) + 0(||uF|li2(p(j))). 
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where x is given by (j3.32|) . Because of (|4.69|) . there holds 
\Up\ - — — 

Fit) 



< 



7s{r) + h 



Fit) 

+0{\\uF\\l2(F(t)-) + WuFWhrgm) 



-fs{r) + h J Fit) 

idn) 



i p. ■iZ- 'jsir) - h 
ls{r) + h 



-ui,x{ls{r) + h)di 



z - 7s(r) - h 



lFit)ls{r)+h ^-^^'"^'^ ■ \^ + 

+0(11^^^1112(^(4)) + 1 1 ^^11 1,2 (9!^) + \\DiuF)\\L^iFit))\\uF\\L2iFit)) 
= 0(^\\uF\\l2(^Fit)) + ll^i^llL2(af7) + \\DiuF)\\L2iFit))\\uF\\L2iFit)))- 

which imphes that 



u^pu%{dnp'^p + d^ipp) 



Fit) 



0[ \\D{uF)\\L2^F(^t))(.\\DiuF)\\ + WuFh^ian)) + WuFWhi^Fit)) 



Combmmg with estimates similar to the previous section for the other terms, we arrive at: 

-HT) 



fj-F I n{s)ds 

h{0) 



<Co{l + T). 



(4.73) 



Introducing the bounds of Lemma 14.21 and Lemma 14.31 into (j3.20p . we show there exists a 
constant C such that: 



n{h)>^, V/iG (0,/iAf), 



so that ()4.73p yields finally : 



ln(/i(T))| <C7o(l+T), VT>0. 



This ends the proof. 
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A Estimates on the test-function in the sUp case 



In this section, we collect estimates on the test-function of Section [3] as defined in (j3.33p . For 
further computations, we define. 



^'(r, z) := $ ( r. 



h + '-fs{r) 

In particular, we note that ^{r, h + 7s(r)) = 1 so that: 



d 
dr 



-^{r,h + -fs{r)) 



2 \ 2 

dr [ y ^ j (r, h + js{r)) + -^^==^d,^{r, h + 7s(r)). 



Thus, on dSh n dQh,i, there holds : 

^2 



n-(ph = ya^^'(r,/i + 75(r)) + ^^^^^a^ ("y^- ) (r,/i + 75(r)). 



yl — r^ = 63 • n . 



Before going to estimates on 99/1, we first obtain a precise decription of the size of ^ and its 
derivatives. Here and in what follows we use a < 6 to denote that a < Cb for a constant C 
independant of relevant parameters. 

Proposition A.l. The following inequalities hold true, for all values of h £ (0, /jm) and 
{r,z) G nh,5 ■ 



\dr,^ir,z) + 



\^{r, z 
\d,^{r,z 

\drr^{r, Z 

-dz'^{r,z 



dzzz'i>{r,z)\ + \drrz'i>{r,Z 



\drzz'^{r,z)\ < — --"2, 



< 1, 



< 



< 



< 



< 



1 



1 
r 

h + ls{r) ' 
1 

\drrr^{r,z)\ < 



\dr^{r.z)\< 
\dzz^{r,z)\< J 



+ 



(/i + 7s(r))2 /i + 7,(r) 



(A.74) 
(A.75) 

(A.76) 

(A.77) 

(A. 78) 

(A.79) 



25 



and also 

\d.nr,z)\ < j-^^, (A.80) 

\drh^{r,z)\ < + ,,, , (A.81) 

(/i + 7s(r)) {h + -ys[r)y 

\d,,h'fir,z)\ + \drrh'i^{r,z)\ + \dh,^{r,z)\ < ——1-—, (A.82) 

\drzh^{r,z)\ < ^ \ ..3 . (A.83) 

Proof. At first, we note that 

\dras\ + \drap\ < r, 

\drrras\ + {drrOsl + \drrrap\ + \drrap\ < 1, 

This yields, by (13:211) - (13:25]) . that we have: 

\drHr, t)\ + \drtt^it, r)\ + r)| < r , V (r, t) G (0, x (0, 1) . 

and 

\dMr,t)\ + \drrrHt,r)\ + |a,rt$(t,r)| + |c>„$(t,r)| < 1. 
We remark also that 

\duMt,r)\ + \dtMt, r)\ < {ap + as) < {h + 7s(r)). 

We introduce these bounds implicitly in the following computations. 
For instance, we have: 

dr^{r, z) = dr^ (r, ^— ^ - —i^ii^l— ( r, 



< 1 

We obtain similarly (IA.75p - (IA.79p after tedious calculations. The only computation we make 
precise is the one concerning 9r-2^- Differentiating the last equality w.r.t. z, we get : 



where 



Thus, 



(/i + 7s(r)) 



r „ , ^ 1 „ r — 7'c.(r) 

-——-d,^{r,z) = .A t'^+ (j,^ )L 

h + -fs{r) [h + jsir]) ih + JsinY 
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where we dominate : 

\r-ls{r)\ ■ 

Combining with the previous estimates on the derivative of we obtain the expected result. 
Concerning derivatives w.r.t. h the computations are similar. First, we note that : 

\d^dhasir)\ + \d^dhapir)\ < 1, Vfc e N. 



Correspondingly, we obtain that 



Then, we obtain, for instance : 

(/i + 7s(r))^ (/i + 75(r))^ 

~ TTT, r~n2 ^i/^^ + TTT, TT^'^"* 

(/i + 7s(r))^ (/i + 75(r))4 

in which, we plug previous estimates on <^ and its derivative. This yields : 



Other estimates are obtained similarly. □ 

The rest of this appendix is devoted to the proofs of lemmas 13.21 and 13.31 They rely on 
the following lemma whose proof can be found in [13] : 

Lemma A. 2. Given {p,q) G (0, oo)^, when h goes to 0, the quantity 

\x\Pdx 
Jo ih + \x\^)i 

behaves like 
p+i 

(1) h—^'i, ifp+l<2q, 

(2) In/i, ifp+l = 2q, 

(3) C, tfp+l>2q, 

with C a constant depending on (p, q) . 

A.l Proof of Lemma 13.21 

First, we recall that, by construction iph is smooth outside ^lh,5- Thus we have get that 

l|VVllLoo(F,An,,,) < 1, V/cGN. 

Hence, we focus on i^h,s and replace iph with (ph in what remains. We keep notations of the 
previous lemma. In particular, we have by definition : 



^ (^-d.-^ir, z)rer + ^dr [r^^^ir, z)] e,^ , in nh,5- 
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Sobolev bounds. First, applying the bounds of the previous lemma, we get: 



1 



< 



1 + 



-5 ^3 



r'^dr 



Then, in cylindrical coordinates, we obtain: 



< 1. 



where we used (|A.77p to get rid of drz"^ terms. Consequently : 



/i + 7^(r)' 



(A.84) 



On the other hand, we have 



in 2 1|2 



^ r2|sin(0)|2 . , ^ 



{h + ls{r)? 
Let now e > small enough. There holds: 



+ 0(1) 



1^) ,^2 ||2 



> J||a.^(r,z)||2,(^^_^)^^|^|^^^ + 0(1) 



> c 



1 rS 



Jo 



\dMr,t)\ 



2 rdr dt 



0(1). 



where c = is a small positive constant. Finally, as dt^ is bounded from below, we get 



in 2 1|2 



> c' 



1 rS 



JO 



rdr dt 



{h + 75 (r)) 



+ 0(1) > c| lnh\. 



Combining the above estimate yields : 



c\lnh\ < ||Z5((^i.)||i2(^(,)) < \\V^F\\l2(^Fit)) < C\lnh\. 
This completes the proof of (j3.36p - (j3.38p . 

Boundary estimates. We proceed with estimate (|3.39p . By standard identities of differ- 
ential geometry (see [3]), we have 

2D{iph)n X n = dn^h xn + {^h- ea) x n. 
Thus, we only need to prove that 

WsdnVh X n + {I3s + l)ifh X ?^llL2{a5hnanh,2i) < C- 
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In the aperture ^h,2S we recall that the normal to dSh directed toward Sh satisfies n 
—rsr + ((1 + h) — z)ez- This yields: 

2^ft X n = ^-(9^^ rer + ^dr p^-] e^^ x {-rcr + ((1 + h) - z)e^) , 

= (r((l + h)- z)d,^ - dr p^] ) ee . 
Note here that, on dSh, we have {1 + h) — z — 1 < r'^ and that 

dr P^] < r = 0(1) . 

Here and in the end of this proof 0(1) means bounded in L'^{dSh H dQh,25) ■ This yields: 

^hxn = ^d.^ee + 0(1) . 
On the other hand, in the aperture ^lh,s we have: 

2dn^h = (rdrzi'^r) - {{1 + h) - z)rd,,^) er 



+ -rdr 



+ ((1 + h)- z)d. 



Introducing the estimates of Proposition lA.il we get in particular: 



9. 



^2 



h + js{r] 



5.^ + 0(1) 



2h 



h + 7s (r) 



5.^ + 0(1). 



Similarly, we have: 

rdrzi'ifr) = r'^drz^ + rdz^ = -rdz^ + r 
using that 



/i + 7s(?') 



+ 2 a,^ = o 



+ 2]dz^ + 0(1) = -rdz"^ + 0(1), 



(A.85) 



rh 



is uniformly bounded in (use case 3 of Lemma lA. 21 with p = 3, g = 4). Then, 



2dn^h = -{rdz-^ + rdzz^)er + 



2/i 



+ 0(1). 



h + 7s (r) 

Recalling that n = —rer + ((1 + h) — z)ez, and still applying (lA.SSp . we obtain: 



Finally, we obtain: 



By our choice of 05, it follows that (see (I3.29P ) 

1 



dzz^ir, h + 7s(r)) + ( — + 2 ) dz^{r, h + 75(0) = 0(1). 



This yields 
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h- derivatives. For dhfh^ we note that, for z = 1, 2, 3 we have: 

Consequently, as integrating with respect to the z-variable reduces to remove one power of 
(h + Jsif)) in the denominator, we get 



h+-y3ir) 



dh^Fds\\L2(^F^nnn,s) ^ ^ + 



{h + -fs{r)) 



< 1, 



and 



II dh^pdzWh^dundn 
Finally, we also compute: for i = 1,2: 



< C|ln/i| . 



|5i9h^h| < |9,/,^'| + |r5,,/,^| + IS./.^'I + |r5,,;,^'| 



< 



1 



(/i + 7s(r))2 
so that for the same reason as above, we have: 

rh+'jsir) 



dh^l{xi,X2,s)ds\ <1, fori =1,2 



and 



+ ls{r))di dh^p{xi,X2,s)ds\\L2(^p^f^Q^^^)<l, fori = 1,2,3. 

J z 

This ends the proof of I^M^-I^M^. 



Pointwise bounds. We easily get that on Q^^s 

\y^h\ < \rdrz^\ + 15,^-1 + |ra,,^| + |a,^| + |ra,,^| < 
from which we conclude: 



1 



h + 7s (r) 



/ {di(p{{xi,X2,s) + dj(fl{xi,X2,s))di 

J z 

Similarly, we get: 



<1, for G {1,2,3}2. 



\dxV^h\ + \d2^^h\ < |r9„,^'| + IS^'I'I + IrS.^.^'l + l^^.^-l + l^^.^-l < — — + 



and we conclude: 

{h + js{r))di J {diip^p{xi,X2,s) + dj(fp{xi,X2,s))ds 

This ends the proof of our lemma. 



< 1. 
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A. 2 Proof of Lemma E3] 

Boundary integrals We focus on the domain fih^s a.nd we replace iph with (^/j in the 
proof because everything is bounded elsewhere. First, it turns out that the first identity is 
independant of qh. Indeed, applying boundary conditions ()3.3ip we obtain : 

{2D{Lph) - qh)n ■ ifhda = / [{2D{ifh)n) x n] ■ {(ph x n)d(J 



1 



Similarly, introducing boundary conditions (|3.3ip and estimate (I3.39p . we compute: 

{2D{^h) - qh)n ■ (es - ^h)<i(T = / [(2L>(^/,)n) x n] ■ ((es - iph) x n)da 



1 

Ts 



1(63 - iph) X n| da + 0(11(63 - iph) x n||i2(a5^nn;, s)) ■ 
dnh,sndSh 



On the boundaries, there holds 

_ Pi{r) ^ T 

which implies that 
and which implies that 



\lnh\< {2D{iph) - qh)n ■ {e-s - iph)da - {D{iph) - qhl)n ■ (fhda < \ln{h)\. 

Jdnh,sr\dSh JdQr\dnh,s 

Construction of q^. Again, we focus on the domain s where iph = 'fh to compute q^. 
The construction is extended to the whole Fh by a standard truncation argument. 

Given /i > 0, let define 



rdrz'^ + 2d,^ + d,,,^{s, z)sds 



Obviously, the map (/i, x) i— )• qh{x) is on {h G (0, /im), x G F^}. and we have Aiph — S/qh 
//i63 where 



\fh\ < Irdrrr'^l + Idrr^^l + 



+ \rdrzzM + \dzzM < ^ 



Then, for arbitrary h E (0, and v G H^{Fh) we introduce x the truncation function as 
defined in (j3.32p and we split : 

/ {AyDh-Vqh)v= / x{^^h-Vqh)v+ / (1 - x)(Av3/, - Vg/,)w = / + //. 
JFh JFh JFh 
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Then, we have: 

\xfh\\v3\ 



V3 



< \\ifT' + 'ys)fh\\mFhnnn,2s)\\^^(^j.^j^^\\L^{Fhnnn,2s) 
where a standard Poincare inequahty yields: 

As for the second integral, we introduce PzR^ for i = 1,2, 3, which satisfy : 

PzR' = J ^^^^ \l - x){^A - ^^qh) for r < 1, z<r + 7,(r) (A.86) 

fh+l 

PzR" = J {l-x)iWh-diqh) otherwise. (A.87) 

Repeating the computations of Ii{t) and l2{t) in Section [3.31 we obtain (the summation over 
index i is implicit) : 

\II\= I {d-iv' + diV^)P,R' + [ v^diP,R'-[ {v\' + v'n^)P,R'da . 

As (A(y3/j — Vqh) is uniformly bounded outside ^h^s, applying a Cauchy Schwarz inequality 
to dominate K2 and combining with Ki concludes the proof : 



/ {IS.iph - ^^Qh) ■ V 



^ \\D{v)\\L^F^,) + \\v\\LHdn)- 



B Estimates on the test function in the mixed case 

In this appendix, we collect estimates on the test function of Section [H As previously, we 
define, 

il,{r,z):=^(r,-—^—], yhe{0,hM), V (r, z) G . 



The only point which diff'ers from the test-function in the slip case is the decription of this 
^. Combining these estimates in the same manner as in the previous appendix, we obtain 
Lemma [4.2[ As the computations are completely similar, all the proofs are left to the reader. 
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Proposition B.l. The following inequalities hold true, for all h G (0, /ijv/) and (r, z) £ ftf^ s '■ 



\drz'if{r,z 



< 1, 



< 



< 



< 



< 



< 



1 



1 

r 

(/i + 75(0)^ 
1 

(/i + 75(r))2 
r 



1 



r 
r 

{h + ls{r)Y 



and also 



\dh^{r, z) 

\drrh'^{r,z) 

\drzh^{r,z) 



< 



< 



< 



1 

1 

{h + ls{r)Y 
r 



+ 



1 



(/i + 75(r))3 (/i + 75(r))2 
1 



(/i + 75(r))3 



B.88) 
B.89) 

B.90) 

B.91) 

B.92) 

B.93) 
B.94) 

B.95) 
B.96) 
B.97) 
B.98) 



We also need a different way to construct the pressure qh in order to prove Lemma l4.3i 
Again, we construct qh only in tlie aperture domain ^h,s-, the extension to being done by 
a standard truncation. So, for /i > we set : 

2qh=rdrz^ + 2dz^ - / d^zz^sds. 

Jo 

In this way, we obtain Acph — S/qh = f^er + /^e^ where the /j's are bounded outside Q,h,5 
and, in fthS^ we have: 



\drn + 



in < \rdrrr'f\ + \drr'f\ + 

in < |r5,,,^| + |5,,^| < 

r 



< 



< 



IrS^^^^^-l + {drrz'ifl + 



< 



(/i + 75(r))2- 

Then, for any v G H^{Q,) such that • n = on dO, and v G 7^ on S/i we have 



(B.99) 
(B.lOO) 

(B.lOl) 



fv^ + I fer ■ V. 
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First, we apply Poincare inequality, as f ^ = on dQ n dQh,s, together with (jB.99p . This 
yields: 



< 



V3 



For the second integral: we write f^Cr = f^ei + /^e2, and split the integral, with x as defined 
in 1^:^ : 



rv^= / X/V + 0(||7;|U2(^^)), 
Fh JFh 

as is uniformly bounded outside Vth 5- Then, we introduce PzP , for i = 1, 2, which satisfies: 



Pzf = / XP 
'0 



Through integration by parts, we get (with the Einstein convention of repeated indices): 



Fh JFhnnh,26 JFhnQh,2s + ^ysK^) 



Pzf{v'n-' + v^n') = I + 11 + III 

idShnUhas 

Thanks to the bounds (IB.lOOD - dB.lOip . we get: 

+ < \\D{v)\\L2^p,^y 

As regards the last integral, we simply write, because ^i^,^ G TZ, 

^ \\Pzr\\L^dSh)M\L^{dSh) ^ lbllL2(5;.) ^ lbllL2(n) 



This ends the proof. 
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